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Abstract 



We present a scenario of neutrino masses and mixing angles. Eacli gen- 
eration includes a sterile right handed neutrino in addition to the usual left 
handed one. We assume a hierarchy in their Dirac masses similar to, but much 
larger than the hierarchies in the quarks and charged leptons. In addition, we 
include a Majorana mass term for the sterile neutrinos only. These assump- 
tions prove sufficient to accomodate scales of mass differences and mixing 
angles consistent with all existing neutrino oscillation data. 
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I. INTRODUCTION 



While there is currently no direct experimental evidence for neutrino masses, there is 
growing indirect evidence in the form of neutrino oscillations, culminating in the recent ob- 
servation of — > Ve at LSND The combined evidence suggests three independent mass 
splittings among the neutrinos participating in the oscillations. If each of these splittings 
is taken seriously, then a fourth neutrino is required to accomodate all the data [0. Sev- 
eral such scenarios have been proposed [Q. In most cases small mass differences and small 
mixing angles are put in by hand, and it seems difficult to explain their origin without fine 
tuning. Furthermore, these scenarios treat generations on an unequal footing, mixing the 
extra sterile state with only the electron neutrino. 

We examine here the viability of one sterile neutrino for each generation. Such models 
have not been considered previously due to the constraint N,j ^ 4 from big bang nucleosyn- 
thesis [Q. It is known, however, that this constraint can be avoided if the tau neutrinos 
have masses in the MeV range and decay rapidly into 0- The three mass splittings 
then suggest a unique natural mass spectrum. When combined with various terrestrial ex- 
perimental data, the solar neutrino deficit implies a neutrino almost degenerate with Ue, 
the atmospheric deficit implies a neutrino almost degenerate with i/^, and the LSND data 
implies the two pairs must be split by at least 0.1 eV. We therefore impose Dirac masses 
with a very large hierarchy and a CKM matrix analogous to the one in the quark sector. We 
then include a Majorana mass matrix on the right handed neutrinos only. The scale of the 
Majorana masses is (9(10^^ eV) and is appropriate for a seesaw mechanism between the 
Grand Unification (GUT) scale and the electroweak scale. Neutrino mixings appropriate 
for the oscillation data will arise from the interplay between the Dirac and Majorana mass 
matrices. 
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II. FERMION MASS HIERARCHIES 



Typically a model with a large mass hierarchy will have mass matrix elements whose 
scales obey 



^ A 3 2 ^ 

e"* 



e 1 



X mo; 



e< 1. 



(1) 



Here rriQ is the scale of the largest mass eigenvalue. If the off diagonal elements are larger 
than these, the lighter masses will recieve large seesaw contributions and the hierarchy 
will be destroyed. Note that these are only upper bounds on the scales of the matrix 
elements. Realistic models typically contain texture zeros [0 or additional powers of e in 
their off diagonal elements. Such suppressions are in fact necessary to make contact with 
the Standard Model, as we will see below. In the Standard Model we have ~ 1/14.3, 
ed ~ 1/5.1, and eg ~ 1/7.6 (see Figure [|). With m^^ in the MeV range, the effects of the 
off diagonal elements of the neutrino Dirac mass matrix will be washed out by the charged 
lepton mixings. 

The Dirac mass term in the neutrino Lagrangian is 



rriij is diagonalized by 
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and the charged leptons are diagonalized by 
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(4) 



where we have suppressed generation indices. Thus for example, the neutrino mass term 
can be written 



(5) 



where U^TTLijU^ = diag(mi, m2, ^3). 

If both the up-hke and down-hke members of an SU (2) muhiplet have mass matrices 
hke (|l|), then the down-hke mass eigenstates wiU be rotated from the weak partners of the 
up-hke mass states by the CKM matrix, 
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where e is the larger of e„ and ea- For example, the weak partner of the electron would be 



V 



(6) 



(7) 



and the CKM hierarchy parameter is eg. We assume the hierarchies of the e and v states 
are aligned with respect to each other. It is possible to consider two hierarchies which are 
related by a generic unitary transformation, but it is hard to imagine a mechanism which 
would generate such hierarchies naturally. In a generic basis, the matrix elements would 
appear to be fine tuned to 0{e^). The alignment of the quark hierarchies is evident in the 
smallness of the off diagonal elements of their CKM matrix. The Cabibbo angle in particular 
is very close to its expected magnitude. Note, however, that the other off diagonal elements 
of are smaller by an additional factor of than would be expected on the basis of the 
quark hierarchies alone. Thus is well described by the Wolfenstein parametrization [Q, 

^ l-eV2 e e'Az^ 
V'^= -e l-eV2 e^A 

e^A{l - z) -e^A 1 

where A ~ 0{1) is real, e ~ is real, and 2: is a complex number with magnitude 0{1). 
A priori there is no way to know whether will follow this pattern or the pattern of 
or some other texture. For definiteness we first use pattern (^) with e ~ eg as an example 
and then discuss other possibilities. 
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III. THE MAJORANA MASS MATRIX 



In order to achieve mass splittings and mixing angles appropriate for the Atmospheric 
and Solar oscillation data, we include Majorana masses of (9(10~^eV) for the right handed 
neutrinos. This is about the right scale to be generated by a seesaw between the GUT 
scale and the electroweak scale, although it is not obvious how sterile particles relate to 
the electroweak scale. It is interesting to note that with this interpretation, the electron 
neutrino will be at the bottom of a two stage seesaw. We write the Majorana mass matrix 
as 



Oil ai2 ai3 



m,- 



(9) 



ai2 022 0-23 
y Oi3 023 O33 

where a priori each a^- is O(10~^eV). 

We are now in a position to write the full 6x6 mass matrix. We write the right handed 
neutrinos in terms of their charge conjugates, vr = s'' . We use (^, @, and with the 
pattern (p|) to get 



= -NCM'N + H.c. 
where C is the charge conjugation matrix, N = (z/g, Se, v^, s^, Ur, St), and 
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The rrii are defined by (|^). 

We take Im(2;) = 0. We will not consider CP violation here. M' is diagonalized by 
M' = OPMPO, where O is orthogonal and P is a diagonal phase matrix. The weak 



eigenstates are written in terms of the mass eigenstates by — Uai Vi — OiaPi vi- We 
have 
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and the masses are 

I — ' «ii ' ^2 - ^ , ^2 + — , m^[l + ^-) , m3(l + I , (13) 

where we have dropped higher order terms. 



IV. EXPERIMENTS 



The probabihty that an initial Va of energy E will oscillate into after a distance L is 



^ ^/3) = <^a/3 - 2 E Re[C/„,C/^C/^,C/;3,] sin^ ^ (14) 



2E ^ ^ 

For the LSND experiment, the relevant terms are 

P^e = e (2sm — + 2sm — sm — ). (16) 

The last term vanishes since its wavelength is too long for the 30m LSND basehne. The 
other terms are effectively equal, and together have the same effect as a two flavor — i/^ 
oscillation with SM"^ = ml and sin^ 2^ = 4e^. We may therefore appeal to published two 
flavor analyses. 
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The LSND experiment was designed to be most sensitive at the mass sphtting preferred 
by the CHDM model of cosmological structure formation 0, ~ 6eV^. The allowed 
regions of sin^ 26 at that depend on how the data is analyzed. The 99% likelihood region 
is 0.002 < sin^ 29 < 0.01, while the 80% confidence level band is 0.0012 < sin^ 29 < 0.005 
. The confidence band uses only the number of events, while the likelihood region uses all 
the information about the events including neutrino energy and distance from source to 
detection point so it is the best way to determine favored regions of and sin^ 2^^ [1T0| . 
The difference is important because 0.002 < sin^ 29 is excluded at 90% confidence level by 
the BNL E776 experiment ||ll|. Thus if we compare similar types of bounds, there is a 
marginally allowed region consistent with the cosmologically preferred 5M'^. 

On the other hand, there is a large region at lower SM"^ allowed by all the data (including 
limits from r-process nucleosynthesis |jl2[)- The 90% likelihood region of LSND, combined 



with the 90% confidence limits from E776 and the Bugey reactor experiment allow 
0.25 < 5M^ < 2.3 eV^ and 0.002 < sin^ 29 < 0.04, giving 0.5 < ms < 1.5 eV and ^ < e < ^. 
The upper end of the range of e is reasonably close to our expectation of e ~ ^ . If sin^ 26'is 
found to be in the lower end of this range, then the mixing of the charged leptons must 
be supressed from its expected hierarchical value for almost any conceivable neutrino mass 
scenario. 

For the atmospheric deficit, the relevant probability is 

i^M. = l-sin^% (17) 

This has the same effect as a two flavor u^—Vs oscillation with = 2m2a22 and sin^ 29 = 1. 
Maximal mixing is allowed by the combined Frejus, NUSEX, 1MB, Kamioka sub-GeV, and 
Kamioka multi-GeV zenith angle dependent data for 4 x 10~^ < < 0.01 eV^ [Q, giving 
0.13 < 022 < lOmeV (milli-eV). There is a small probability (~ 2%) of oscillation into u^. 
Matter effects are insignificant. 

The vacuum disappearance probability for Ug is 

p;- = 1 _ 2e2 sin2 ^ _ 2e^ sin^ (18) 
7 
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FIG. 1. Fermion mass hierarchies. The up and down quark and charged lepton hierarchies are 
shown. We include the proposed Dirac neutrino masses. The bars indicate the range of masses 
consistent with experiments. 

This is too small an effect for the solar neutrino deficit, so we assume a. Pe — small 
angle MSW mechanism. The allowed parameters are 0.003 < sin^ 29 < 0.012 and 4 < 
< llmeV^. We have SM'^ = aj^ and sin 6' = giving 2 < an < 3.3 meV and 
0.05 < nil < 0.2 eV . Note the similarity in the scales of an and a22- 

The Dirac neutrino masses are plotted in Figure |l]. We note in passing the possibility 
of extending the hierarchies to a fourth generation. Three of the particles would have 
suggestively similar masses. It turns out that to be consistent with weak neutral current 
data [|r^, the fourth up- type quark would have to have a similar mass to the other three 



8 



particles, in conflict with the existing hierarchy. 

The probabilities for appearance experiments are 

p^^ = eM2(- sin^ ^ + sin^ ^ + sin^ ^ + sin^ ^ + sin^ ^ - sin^ ^) 
P^^ = eM2( - 2{z - l)(sin2 ^ + sin^ ^) - 2z{z - l){sm'' ^ + sin^ ^) (19) 
+z{sm'^ ^ + sin^ ^ + sin^ ^ + sin^ ^) - sin^ ^ + z"^ sin^ ^ 
Each expression contains terms where A^- will be large for any conceivable experiment. For 
those terms, sin^ ^ will average to 1/2 over the finite E and L resolution of an experiment. 
Thus we may estimate the probabilities as ~ and Pgr ~ e^- Currently the best limits 
for large are P^r ^ 0.002 and Pgr ^ 0.073 from the E531 experiment at Fermilab [[16 



Pfj^r gets its scale from Vj^ ~ e^. Thus, while the parametrization @ is viable, maximal 
hierarchical mixing with Vj^ ~ e is ruled out for any scenario with 6M^^ > 10 eV^. Maximal 
hierarchical mixing is not constrained in the e-r channel. 

Upcoming experiments may be able to distinguish this scenario from the other possibil- 
ities in the next few years. The prediction that the atmospheric deficit is caused by yU — s 
oscillations with maximal mixing is unique to this scenario. The mixing angle could be 
pinned down with further atmospheric neutrino data. Chooz and San Onofre can 



eliminate the possibility that the atmospheric deficit is fi-e by directly measuring z/g disap- 
pearance probabilities. And ICARUS [||] and MINOS @ might be sufficient to rule out 



yU-r. This would leave fi-s as the only alternative. An MeV u^- is then almost inevitable to 
save BBN since there would be at least four active neutrino flavors at the time of nucleosyn- 
thesis. Observations of the solar neutrino spectrum can firmly establish the mass splitting 
and mixing angle for the neutrinos responsible, and KARMEN |^ can confirm the LSND 



result, which would eliminate possibility for the solar deficit since the mass splitting 

would be too large. An observation of fi-r oscillations at CHORUS [^], NOMAD ||23|, or 



COSMOS 1^ would then firmly establish fi-s for the solar channel. While CHORUS and 
NOMAD themselves have a chance of observing fi-r, COSMOS is very likely to observe this 
channel, but very unlikely to see e-r. It is clear that the next few years will be very exciting 
for neutrino physics. 



V. CONCLUSIONS 



We have shown how neutrino masses appropriate for the various oscillation data can be 
fit into a hierarchical mass scenario analogous to the hierarchies in the quark and charged 
lepton sectors. Small mass splittings and small mixing angles result from the interplay of 
Dirac and Majorana mass terms. The scenario satisfies all experimental and astrophysical 
constraints. It is unique among proposed solutions in that the atmospheric oscillations are 
^-s with maximal mixing, a prediction which could be tested experimentally in upcoming 
experiments. 
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